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Abstract. Classical uniaxially anisotropic Heisenberg and XY antiferromagnets in a field along the easy 
axis on a square lattice are analysed, applying ground state considerations and Monte Carlo techniques. 
The models are known to display antiferromagnetic and spin-flop phases. In the Heisenberg case, a single- 
ion anisotropy is added to the XXZ antiferromagnet, enhancing or competing with the uniaxial exchange 
anisotropy. Its effect on the stability of non-coUinear structures of biconical type is studied. In the case 
of the anisotropic XY antiferromagnet, the transition region between the antiferromagnetic and spin-flop 
phases is found to be dominated by degenerate bidirectional fluctuations. The phase diagram is observed 
to resemble closely that of the XXZ antiferromagnet without single-ion anisotropy. 
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1 Introduction 

Two-dimensional uniaxially anisotropic Heisenberg anti- 
ferromagnets in a field along the easy axis have received a 
renewed interest in recent years, experimentally as well as 
theoretically. On the experimental side, especially layered 
cuprates exhibit interesting properties due to an interplay 
of spin and charge [Il[ll[31lll[51ini[71[H] . Intriguing phase dia- 
grams have been obtained for other quasi two-dimensional 
antiferromagnets as well, showing, typically, multicritical 
behaviour [gifigilTTIfT^ . 

On the theoretical side, recent studies [l3 l fT4 l [T5 p i6 ( fT7| . 
[T51ll9'"20' on the square lattice XXZ Heisenberg antiferro- 
magnet have substantially extended previous analyses [21|, 
[22 ll23[l24| on this prototypical model. The XXZ model is 
described by the Hamiltonian 

Hxxz = J E [ A{S^S^ + SfS]) + S^S^] - S! 

id i 

(1) 

where the sum runs over all pairs of neighbouring sites, 
i and j, of the lattice, J(> 0) is the coupling constant, 
and A is the exchange anisotropy parameter, with A=0 
corresponding to the Ising limit and A~ 1 to the isotropic 
Heisenberg case. H is the external field along the easy axis, 
the z-axis. S", a=x, y and z, are the three components of 
classical or quantum spins. 

For many years |21| . the model is known to display, 
at low temperatures and small fields, a long-range or- 
dered antiferromagnetic (AF) phase, and, at higher fields 
and low temperatures, a spin-flop (SF) phase with alge- 
braically decaying correlations. However, the multicritical 
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Fig. 1. Antiferromagnetic (a), spin-flop (b,d), and biconical 
(c) as well as bidirectional (e) structures which may occur 
as ground states in the XXZ antiferromagnet and variants, 
(a,b,c), or in the anisotropic XY antiferromagnet (a,d,e). The 
tilt angles are defined with respect to the easy axis, being the 
z-axis in the XXZ and the x-axis in the XY case. 



point, where the AF, SF and paramagnetic phases meet, 
had been subject to controversial suggestions. 

In early renormalisation group calculations [25 p 26 ( l27] 
the XXZ Heisenberg model and its extensions to n-compo- 
nent antiferromagnets with uniaxial anisotropy have been 
investigated. Various possible multicritical scenarios have 
been proposed, depending on the number of spin compo- 
nents, n, and the dimension of the lattice, d. The possible 
scenarios include a bicritical point of 0(n) symmetry, a 
tetracritical point, and a critical end point. Actually, an 
e-expansion to low order, favours, for n = 3, i.e. for the 
XXZ antiferromagnet, the bicritical point. Certainly, the 
bicritical point can not be realized in two dimensions at 
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a non-zero temperature, T > 0, because it would violate 
the rigorously proven, well-known theorem of Mermin and 
Wagner j8j. 

Recent Monte Carlo studies and ground-state consid- 
erations provide strong evidence for the multicritical point 
being a 'hidden tetracritical point' at T=0 in the classical 
XXZ model P^IT^ITT] . Indeed a narrow phase governed 
by 'biconical' [27] (BC) fluctuations separates the AF and 
SF phases at low temperatures. At that zero temperature 
special point, AF, SF, and BC structures have the same 
energy, leading to a high degeneracy [17'. Actually, the 
importance of the non-coUinear BC structures, see Fig. 1, 
for the ground state and the phase diagram of the classi- 
cal XXZ antiferromagnet had been overlooked in previous 
work. On the other hand, already a few decades ago, it had 
been noticed that such BC structures may be stabilised by 
adding to the XXZ model further anisotropy terms, like 
cubic terms, or longer-range interactions PW5D1I5T] . 

Extending significantly our very recent short commu- 
nication [20], the aim of the present article is twofold: 
Firstly, we shall study the effect of adding a single-ion an- 
isotropy to the XXZ model, which may either enhance the 
uniaxial exchange anisotropy. A, or it may compete with 
it by being a planar anisotropy, depending on the sign of 
the coupling strength of the single-ion anisotropy. In both 
cases, the special point of high degeneracy (the hidden 
tetracritical point) does not survive. Resulting phase dia- 
grams will be determined using Monte Carlo techniques. 
Secondly, we shall consider the two-component, n — 2, 
variant of the XXZ antiferromagnet. In that case, in prin- 
ciple, a bicritical point, of 0{2) symmetry, would be al- 
lowed to occur at T > 0, being, in two dimensions, in 
the Kosterlitz-Thouless universality class [32]. Note that 
now BC structures are replaced by 'bidirectional' (BD) 
structures, see Fig. 1. Again, in addition to ground state 
considerations, the phase diagram will be determined us- 
ing Monte Carlo simulations. 

The paper is organised as follows: In the next section, 
the models will be introduced and ground state properties 
will be discussed, emphasising the role of non-coUinear, 
BC and BD, structures. Phase diagrams and critical prop- 
erties, as obtained from large-scale simulations, will be 
presented in the then following section. A short summary 
concludes the article. 



2 Models and ground state properties 

In order to study the impact of biconical or bidirectional 
structures on phase diagrams of uniaxially anisotropic an- 
tiferromagnets on a square lattice, we consider two dif- 
ferent classical models with spins of length one. Firstly, a 
single-ion anisotropy is added to the XXZ model, eq. (1), 
so that the Hamiltonian reads 

nD^nxxz + Dj2is-f (2) 

where the single-ion anisotropy may, depending on the 
sign of D, enhance the uniaxial exchange anisotropy A 
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Fig. 2. Ground state phase diagram for A — 0.8. The filled 
circle denotes the highly degenerate point at D = and H = 
2.4J. 

(0 < Z\ < 1), when D < 0, or it may introduce a compet- 
ing planar anisotropy, D > 0. Secondly, the anisotropic 
XY antiferromagnet is studied, described by the Hamilto- 
nian 

nxY = JY.[ s^s^ + Asys^^ ] - ^ 5f (3) 

being the two-component, n = 2, variant of the XXZ an- 
tiferromagnet, where the s-axis is the easy axis. 

The ground state configurations, at T = 0, are fixed 
by the spin orientations on the two sublattices, A and B, 
formed by neighbouring sites of the square lattice. The 
configurations may be determined in a straightforward 
way 29,30,33 . 

For the XXZ model with single-ion anisotropy, eq.(2), 
the xy-components of the spins order antiferromagneti- 
cally, having rotational symmetry. The orientations of the 
spins on the two sublattices are then given by their tilt 
angles, 0a and 6*3, with respect to the easy axis, the 
z-axis, see Fig. 1. From minimisation of the energy, the 
actual values of the tilt angles in the ground state config- 
urations follow as a function of the anisotropy parameters 
A and D as well as the field H . For calculational conve- 
nience, one may substitute Oa and Ob by combinations 
of the z-components of the sublattice magnetisations (per 
sublattice site), Sa = cos 0a and — cos 0b ■ 
For negative couplings D < 0, only AF, SF, and and ferro- 
magnetic (F) ground states occur, see Fig. 2 (setting A= 
0.8, as before [l4l[T5l[T7ll20l[2T] V As depicted in the figure, 
AF ground states are stable for low fields, H < Hd, while 
F ground states are obtained for high fields, H > Hc2- 
The two critical fields are given by, for D > —2AJ, 

H,i = 2v/(2J)2 - i2JA + D)^ (4) 
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Fig. 3. Binder cumulants U^l^ of the longitudinal staggered 
magnetisations of the XXZ and anisotropic XY antiferromag- 
nets at the highly degenerate point T = Q and Hc\ ■ 



and 



Hc2 = 4J(1 + A) + 2D 



(5) 



For intermediate fields H, Hci < H < Hc2, the ground 
states comprise SF configurations, but they are squeezed 
out for a more negative single-ion term, D < ~2AJ, fos- 
tering the uniaxial alignment of the spins, see Fig. 2. 

Note that for D < 0, in the limit A = 1, when the uni- 
axial anisotropy of the Hamiltonian is solely due to the 
single~ion term, no biconical structures exist as ground 
states. In contrast, when the uniaxial anisotropy is solely 
due to the exchange anisotropy, i.e. in the XXZ antifer- 
romagnet, D — 0, BC structures are ground states at the 
critical field H^i = 4JV1 - Z\2, see eq. (4) ^7M\. At this 
highly degenerate point the BC structures take on tilt an- 
gles interrelated by [T7] 



Of 



Vl - A"^ - cosOa 
1 - Vl - A"^ cos Oa 



(6) 



interpolating continuously between the AF and SF config- 
urations, with the tilt angle 0a ranging from to tt. One 
may calculate ground state values of various quantities at 
the degenerate point, assuming that each degenerate con- 
figuration with the interrelated tilt angles occurs with the 
same probability and taking into account the rotational 
invariance of the spin components in the xy— plane. 
As an example we mention the Binder cumulant |34j 



U^, = l-< {Ml,f > /(3 < {M^,r >) 



(7) 



where the longitudinal staggered magnetisation is de- 
fined by Af|t = {S% - S'|j)/2. The cumulant at the highly 



degenerate point {T 
shown in Fig. 3. For 



- 0,H — Hci) as a function of A is 
instance, in the much studied case 



A = 0.8, one obtains = 0.3777... [33]. In simulating 
small systems at low temperatures on approach to the de- 
generate point for various values of A, we confirmed the 
(numerically) exact result for the cumulant. 

An interesting quantity for characterising the BC struc- 
tures is the probability of finding the tilt angle 6*, p{0). 
At the degenerate point, one gets [33] 



piO) 



a sin O 

1 — a cos 



In 



1 + a 
I — a 



1 -1 



(8) 



where a — \/l — A"^ /2. Obviously, there is no full 0(3) 
symmetry, as one may expect for a bicritical point of 
Heisenberg type. One may reproduce that form of p{0) 
by simulating, again, small systems at low temperatures 
on approach to the highly degenerate point. 

We now consider the case of a single-ion anisotropy 
favouring a planar anisotropy, D > 0, competing with the 
uniaxial exchange anisotropy A. Now, biconical structures 
may become ground states in a finite, non-zero range of 
fields, Hcia < H < Hcu, as depicted in Fig. 2. The critical 
field between the AF and BC ground states, Hda, is given 

by 



(9) 



and the upper critical field, -ffcib, separating the BC and 
SF ground states, is given by 



Hcib = A[AJ^- {2 J A + Df ] /H, 



(10) 



At = 2 — 2A, the two critical fields approach zero, and, 
at larger planar anisotropics, there exists no AF ground 
state. 

Of course, the degeneracy in the BC structures occur- 
ring at the special point {T = Q, H = Hd) of the XXZ 
model is now lifted, with the tilt angles, 0a, 0b, chang- 
ing now continuously with the field i?, starting with the 
AF structure at Hda and ending with the SF structure at 
Hcib- The magnetisations on the A and B sublattices in 
this biconical region are interrelated by [33j 



S A S f 



(2JZ\)2 

i2j + Dy 



(1 + s\si 



(11) 



The expression transforms into eq. (6) for vanishing single- 
ion anisotropy, D = 0. 

The ground state properties of the anisotropic XY an- 
tiferromagnet, eq.(3), are closely related to those of the 
XXZ model, eq. (1). There is also a high degeneracy at 
the critical field separating the AF and SF ground states, 
Hci/ J = 4\/l — A^, in this case due to bidirectional struc- 
tures. As depicted in Fig. 1, the BD configurations are 
characterised by tilt angles, 0a and 0b, which are again 
defined with respect to the easy axis, being now the x- 
axis. Note that here the tilt angles may vary from to 27r. 
At the highly degenerate point, 0a and 0b are interre- 
lated analogously to eq. (6) . 
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At T = and H = Hd, various quantities of interest 
may be calculated, assuming again that each configura- 
tion with interrelated tilt angles, being an AF, a BD or a 
SF state, occurs with the same probability. In contrast to 
the XXZ case, there is now no planar rotational invariance 
to be taken into account. For example, the dependence of 
the Binder cumulant of the longitudinal staggered mag- 
netisation, C/Jj, on the exchange anisotropy A is included 
in Fig. 3. The probability ^(6*) for encountering the tilt 
angle is found to be 



p{0) = 



1 



1 — a cos 27r 



(12) 



where a = Vl — ^^/2. Analogously to the XXZ case, 
there is no full 0(2) symmetry. Again, these ground state 
results have been checked in simulations, as for the XXZ 
case. 



3 Phase diagrams 

To study the effect of the presence of BC or BD structures 
on the phase diagrams of uniaxially anisotropic antiferro- 
magnets on a square lattice, we shall consider the XXZ 
Heisenberg antiferromagnet with an additional single-ion 
anisotropy, eq. (2), as well as the anisotropic XY antifer- 
romagnet, eq. (3). In all cases the exchange anisotropy 
anisotropy is set A = 0.8, as before. 

Large-scale Monte Carlo simulations have been per- 
formed, studying lattices with up to — 240^ spins, 
with runs of typically up to, for larger lattices, 10* Monte 
Carlo steps per spin, averaging over several realizations to 
estimate standard deviations for the computed quantities. 
These quantities include the specific heat, sublattice, stag- 
gered, and total magnetisations, longitudinal (relative to 
the direction of the applied field) and transverse staggered 
susceptibilities. Binder cumulants and related histograms. 
To monitor BC and BD fluctuations and structures, we 
recorded probability functions of the tilt angles, such as 
the probability ^2(6*^! 0b) for finding the two angles, 0a 
and (9b, at neighbouring sites and the probability ^(6*) 
for encountering the tilt angle 0. The probabilities may 
be defined 'locally' by taking into account the orienta- 
tions of individual spins, or 'globally' by computing av- 
erage spin orientations on the two sublattices. Obviously, 
the two definitions coincide at zero temperature. 

To determine transition temperatures, the finite-size 
dependence of various quantities has been recorded, with 
reasonable extrapolations to the thermodynamic limit (the 
exact finite-size dependence is not always known) . The es- 
timates we obtained agreed within the error bars shown 
in the phase diagrams. 



3.1 XXZ antiferromagnets with single-ion anisotropy 

Let us consider first the case of a negative single-ion an- 
isotropy, D < 0, enhancing the exchange anisotropy. In 
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Fig. 4. Phase diagram of the XXZ antiferromagnet with a 
single-ion anisotropy, A = 0.8 and D/J — —0.2. 



this case, there are no ground states of BC type. In Fig. 4, 
a typical phase diagram is depicted, where D/J= —0.2. 
At low temperatures, we observe a transition of first or- 
der separating the AF and SF phases. Evidence for that 
kind of phase transition has been presented before |20j . 
For instance, the maximum of the longitudinal staggered 
susceptibility is seen to grow with system size, L, like L", 
and a approaches very closely 2 already for quite small 
sizes, L > 10. That behaviour is characteristic for a first- 
order transition with a rather small correlation length at 
the transition. 

Further evidence for a first-order transition may be in- 
ferred from the probability p{0) for finding the tilt angle 
0. Close to the transition, p{0) shows more and more pro- 
nounced local maxima simultaneously at the values of 
characterising the AF phase as well as the SF phase, when 
increasing the system size. Note that at least for small sys- 
tem sizes, biconical fluctuations are also observed in the 
transition region between the AF and SF phases, but the 
relevant effect seems to be the coexistence phenomenon. 

To identify the nature of the triple point, where the 
AF, SF, and paramagnetic phases meet, a very fine reso- 
lution, in temperature and field, is needed near that point. 
This aspect, requiring, presumably, huge computational 
efforts, is beyond the scope of the present study. 
When applying a planar single-ion anisotropy, D > 0, 
competing with the uniaxial exchange anisotropy A, a 
very different phase diagram results. An example is shown 
in Fig. 5 for £)/ J = 0.2. In this case, biconical structures 
occur as ground states in a non-zero range of fields, bor- 
dered by Hcia and Hdb, see eqs. (9) and (10). At T > 0, 
they are expected to give rise to an ordered BC phase, in 
which the ordered AF and SF phases coexist [251130] . Ac- 
tually, here in two dimensions, the algebraic order of the 
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Fig. 5. Phase diagram of the XXZ antiferromagnet with a 
competing single-ion anisotropy, A = 0.8 and D/J = 0.2. 




10 20 40 80 



160 



Fig. 6. Critical exponent for staggered susceptibility x for the 
XXZ antiferromagnet with a competing single-ion anisotropy, 
at fcsT/J = 0.2, see Fig. 5. 



SF phase, is found to vanish at the boundary of the AF 
and BC phases, at Hda, while M|j, the order parameter 
of the AF phase, vanishes at the higher critical field Hcu, 
separating the BC and SF phases, compare to Fig. 5. 

Based on renormalisation group calculations [35,!36ll37l 
[38] , the transition between the BC and SF phases may be 
argued to be in the Ising universality class, while the tran- 
sition between the BC and AF phases is expected to be in 




Fig. 7. Histograms for the probability of the tilt angle p{0) 
for the XXZ antiferromagnet with a competing single-ion an- 
isotropy, D/J = 0.2, at ksT / J = 0.2, at the fields given in 
the inset. Lattices with SO'^ spins are simulated. Note that the 
area under the curve has been normalised to be one. 



the XY universality class, being the Kosterlitz-Thouless 
universality class [32] in two dimensions. 

This description is in accordance with our simulational 
data. For instance, we monitored the size-dependence of 
the maximum of the longitudinal staggered susceptibil- 
ity, Xmax{L), being located close to i?ci&, see Fig. 6 for 
ksT/ J = 0.2. From the doubly logarithmic plot shown 
in that figure, one observes that the effective exponent 
a, defined by Xmax oc seems to approach, indeed, the 
asymptotic Ising value of 7/4 for rather large system sizes, 
L > 120. Thence, significant corrections to scaling play an 
important role. In turn, at the boundary line between the 
BC and AF phases the algebraic order in the transverse 
staggered magnetisation, which holds in the BC phase, 
gets lost. The finite-size dependence of that magnetisa- 
tion, for i > 40, agrees with the transition belonging to 
the Kosterlitz-Thouless universality class, where the or- 
der parameter vanishes at the transition in the form of a 
power-law with an exponent 77 = 1/4. 
In the BC phase the dominant interrelated tilt angles are 
changing continuously, at fixed low temperature, with the 
field. In fact, this behaviour is displayed by the probabil- 
ity function p{0), as illustrated in Fig. 7 for the global 
(sublattice) spin orientations. By increasing the field, at 
ksT/J — 0.2, the peak positions correspond first to the 
AF structure, shifting gradually towards each other, re- 
flecting BC structures, and finally collapsing in one peak 
characterising the SF phase. 

As seen in Fig. 5, the extent of the BC phase shrinks 
with increasing temperature. Eventually, the BC phase 
may terminate at a tetracritical point [5D1I551I55 11371155] . 
where the AF, SF, BC, and paramagnetic phases meet. 
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Fig. 8. Phase diagram of the anisotropic XY model with A — 
0.8. 

Because the phase boundaries seem to meet there with 
common tangents, we give here only a rough estimate for 
the case depicted in Fig. 5 (D/J = 0.2 and A = 0.8), 
ksTtetral J = 0.35 ± 0.05. A more precise location and an 
analysis of its critical properties seem to require enormous 
computational efforts. 

3.2 Anisotropic XY antiferromagnet 

Finally, let us consider the anisotropic XY model, setting 
the exchange anisotropy Z\= 0.8. Its phase diagram is de- 
picted in Fig. 8. 

The topology of the phase diagram looks like in the 
XXZ case [14,15,17 . The AF and SF boundary hues ap- 
proach each other very closely near the maximum of the 
AF phase boundary in the (T, 7?)-plane. Accordingly, at 
low temperatures, there seems to be either a direct tran- 
sition between the AF and SF phases, or two separate 
transitions with an extremely narrow intervening phase 
may occur. At zero temperature and H — Hci{— 2. 4 J), 
the highly degenerate ground state comprises SF, AF, and 
bidirectional structures. 

Away from that intriguing transition region, see Fig. 
8, one expects the transitions between the paramagnetic 
and the AF as well as the SF phases to be in the Ising 
universality class, because in the SF phase of the XY anti- 
ferromagnet, there is just one ordering component, the y- 
component. This consideration is confirmed by the Monte 
Carlo data for the specific heat (where the peak at the AF 
phase boundary gets rather weak on approach to the tran- 
sition region) and for the staggered susceptibilities. The 
quantities exhibit critical behaviour of Ising-type, as fol- 
lows from the corresponding effective exponents describing 
size dependences of the various peak heights. 




Fig. 9. Probability P2(^A,6's) for the anisotropic XY antifer- 
romagnet with A= 0.8 for a system with 100 x 100 lattice sites 
near the transition region between the AF and SF phases at 
ksT/J = 0.4 and H/J = (a) 2.40, (b) 2.4136, (c) 2.4138, and 
(d) 2.43. P2{0a,0b) is proportional to the grayscale. The su- 
perimposed solid line depicts the relation between the two tilt 
angles in the ground state. 



In the transition region of the AF and SF phases, BD 
fluctuations dominate, as one may conveniently infer from 
the (local) probability distribution ^2(6*71, 6*5) for finding 
the two tilt angles at neighbouring sites, i.e. for the two 
different sublattices. Typical results are depicted in Fig. 9, 
showing the behaviour oip2 in a grayscale representation, 
at fixed temperature, ksT/J — 0.4 and different fields, 
somewhat below and above the, possibly, two closeby tran- 
sitions as well as in their immediate vicinity. Indeed, at the 
low field, see Fig. 9(a), one finds a distribution correspond- 
ing to the AF phase, with peaks at Oa/b = 0, tt. At the 
high field. Fig. 9(d), there is a distribution typical for the 
SF phase, with peaks at the spin-flop tilt angle. Note that 
the local maxima in p2 follow closely the line describing 
the dependence of the two tilt angles in the degenerate 
ground state, eq. (6), albeit the probability of these BD 
structures is low compared to those of the AF and SF, 
resp., structures, see Figs. 9(a) and 9(d). In other words, 
thermal fluctuations driving the system away from the AF 
or SF structures are rather weak and, predominantly, of 
BD type. 

In the immediate vicinity of the transitions, see Figs. 
9(b) and 9(c), the BD fluctuations and structures clearly 
dominate. Now, all those degenerate bidirectional struc- 
tures occur simultaneously with (almost) equal probabil- 
ity, i.e. along the line of maxima p2 is (almost) constant. 
Actually, it is also interesting to monitor the time evolu- 
tion of Monte Carlo configurations and of the probability 
p{0) in those configurations. One notices that at a given 
time of the simulation a concrete combination of interre- 
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lated tilt angles prevails. As Monte Carlo time evolves, 
other combinations in accordance with the ground state 
degeneracy prevail, leading to the behaviour depicted in 
Fig. 9. Obviously, this is is marked contrast to the sit- 
uation in the ordered biconical phase at fixed field and 
temperature. There, after equilibration, just one combi- 
nation of tilt angles seems to dominate during the entire 
simulation, see Fig. 7. 

Note that in the transition region at higher tempera- 
tures, analysing systems of fixed size, say, L = 100 (see 
Fig. 9), the various BD structures tend to occur simul- 
taneously in a single MC configuration. Thence, for each 
configuration, p{0) has a shape quite similar to that in 
the degenerate ground state. This observation may be ex- 
plained by a smaller correlation length. Typically, when 
increasing the system size, the region of (almost) constant 
values of p2 along the line of local maxima shrinks some- 
what. Certainly, further systematic studies of the time 
scales and temperature as well as size dependences related 
to the BD structures would be desirable. 

Critical phenomena in the transition region between 
the AF and SF phases have been studied also by analysing 
effective exponents. We did that for the staggered suscep- 
tibilities and the specific heat. Results are compatible with 
Ising-type criticality, but rather large corrections to scal- 
ing had to be presumed. For example, at H/ J = 2.44 and 
0.54 < ksT/J < 0.57, see Fig. 8, the effective critical ex- 
ponents for describing the size-dependences of the peak 
height for the staggered susceptibilities are about 1.8 to 
1.85, largely independent of system size. The supposedly 
rather strong corrections to scaling may be due to very 
large correlation lengths in that region, and the asymp- 
totics may be reached only for very large systems. 

In any event, the simulational data seem to suggest for 
the anisotropic XY antiferromagnet the existence of an 
extremely narrow, disordered phase, intervening between 
the AF and SF phases, hke in the XXZ case [l4llT5] . That 
intermediate phase is dominated by all the, in the ground 
state completely degenerate, bidirectional fluctuations. 

4 Summary 

In this article, we studied two variants of the XXZ an- 
tiferromagnet in a field along the easy axis on a square 
lattice, by, firstly, adding a single-ion anisotropy, and by, 
secondly, reducing the number of spin components to two, 
yielding the anisotropic XY antiferromagnet. Large-scale 
Monte Carlo simulations have been performed, augmented 
by ground state calculations. 

Adding a single-ion anisotropy has a drastic impact 
both on ground state properties and the phase diagram. 
Biconical structures, leading to a highly degenerate ground 
state in the XXZ antiferromagnet, are either suppressed, 
when the single-ion anisotropy fosters the uniaxial ex- 
change anisotropy, or their degeneracy will be lifted by 
stabilising them successively with changing field, when the 
single-ion anisotropy introduces a planar anisotropy. In 
the former case, we observe, at low temperatures a direct 
first-order transition between the AF to SF phases, while 



in the other case, an ordered biconical phase emerges, sep- 
arating the AF and SF phases. The situation in the XXZ 
case without single-ion anisotropy, where a narrow disor- 
dered phase with biconical fiuctuations separates the AF 
and SF phases, interpolates in between these two scenarios 
with single-ion anisotropies of different sign. 

When the uniaxiality of the antiferromagnet is solely 
due to a single-ion anisotropy, i.e. when the exchange 
couplings are isotropic, no biconical structures occur as 
ground states. 

Ground state properties and the phase diagram of the 
anisotropic XY antiferromagnet are observed to resemble 
rather closely those of the XXZ antiferromagnet. There is 
a highly degenerate ground state, at which non-coUinear 
structures of bidirectional type become stable. These de- 
generate bidirectional structures prevail at low temper- 
atures in the transition region between the AF and SF 
phases, leading, presumably, to a very narrow intervening 
disordered phase. 

We conclude that ground state properties and phase 
diagrams of classical uniaxially anisotropic antiferromag- 
nets in two dimensions depend crucially on the form of the 
anisotropy terms, supporting or suppressing non-coUinear 
structures of biconical or bidirectional type. 
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and discussions. We gratefully acknowledge financial sup- 
port by the Deutsche Forschungsgemeinschaft under grant 
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